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In mean-field descriptions of nodal d-wave superconductors, generic edges exhibit dispersionless 
Majorana fermion bands at zero-energy. These states give rise to an extensive ground-state degen- 
eracy, and are protected by time-reversal (TR) symmetry. We argue that the infinite density of 
states of these flat-bands make them inherently unstable to interactions, and show that repulsive 
interactions lead to edge FM which splits the flat bands. This edge FM offers an explanation for the 
observation of splitting of zero-bias peaks in edge tunneling in High-Tc cuprate superconductors. 
We argue that this mechanism for splitting is more likely than previously proposed scenarios, and 
describe its experimental consequences. 



The discovery of topological insulators [1] has led to a 
re-examination of the role of symmetries for protecting 
surface states. Since known topological invariants are de- 
fined only in the presence of an energy gap, this effort has 
focused almost exclusively on gapped insulators and su- 
perconductors. However, well-defined topologically pro- 
tected surface states can emerge at the boundary of (non- 
interacting) gapless systems, so long as translational in- 
variance along the boundary is preserved. For example, 
in the absence of interactions, clean systems with bulk 
Dirac nodes, such as graphene[2], Weyl-semimetals[3, 4], 
and nodal superconductors [5, 6] are all expected to ex- 
hibit dispersionless bands that are spatially confined to 
the edge. In special cases flat-bands can persist in gapless 
superconductors even in the presence of disorder [7]. 

These flat edge-bands exist in regions of linear size A 
in the boundary-Brillouin zone, and terminate at bulk 
gapless nodes where the edge states delocalize from the 
boundary. In superconducting cases, these edge bands 
are pinned to the Fermi-energy as neutral Majorana 
fermions, giving rise to large ground-state degeneracy: 
D 2( A / 27r ) , where d is the spatial dimension of the 
bulk. This corresponds to an extensive ground-state en- 
tropy: So ~ 1 in violation of the third law of 
thermodynamics, and also to an infinite density of states 
at zero-energy. Therefore even arbitrarily weak interac- 
tions act as a singular perturbation, and play an essential 
role in determining true edge-state structure. 

In this paper, we examine the effects of interactions 
on flat edge-bands, focusing specifically on the case of 
2D superconductors with nodal d x 2_ y 2 pairing symme- 
try, relevant to the cuprate family of high-temperature 
superconductors (SC). The flatness of the edge-states is 
protected not only by translation symmetry along the 
edge, but also time-reversal symmetry (TRS). Therefore, 
TRS breaking order is required lift the edge degeneracy. 

Evidence for such TRS breaking was found in normal- 
metal/YBCO tunnel junction experiments. These show a 
large zero-bias peak at intermediate temperatures[8-10], 
corresponding to many low-lying states. The peak sub- 
sequently splits into two as the sample is cooled[ll, 12]. 




FIG. 1. a) Momentum resolved spectrum for the edge for non- 
interacting d-wave nodal superconductor for t — 240meV, 
Ao = lOmeV, and U — (black-line, non-interacting) and 
U = 2eV (red curves, mean- field). Solid lines indicate edge 
states and gray, shaded regions represent the bulk continuum. 
In the presence of interactions the edge becomes ferromag- 
netic, splitting the flat band of edge-states. The small split- 
ting at k — is a finite size artifact. Inset shows the 2D Bril- 
louin zone, rotated coordinate system, and projection onto 
the edge, b) Edge tunneling conductance for 70 = 0.1^-, for 
U = (black) and U = 2eV (red). 

Previously proposed explanations[13, 14] of this phenom- 
ena, were based on the point of view that the edge- 
scattering is pair-breaking and suppresses the d x i_ y i pa- 
rameter near the edge. In principle, this then allows 
a different, sub-dominant pairing symmetry to develop 
near the edge. In particular, s-wave pairing with a rel- 
ative ^ phase to the bulk order parameter (d + is) was 
suggested to develop near the edge[13, 14]. However, 
the relevance of this scenario to the high-Tc cuprates is 
questionable, since superconductivity in these materials 
is widely believed to arise from strong electron-electron 
repulsion[15-17], which disfavors s-wave pairing[18]. 

Viewing the zero-bias peak, instead, from the perspec- 
tive of topological flat bands naturally suggests a dif- 
ferent scenario. Namely, since the flat edge-bands are 
spin-degenerate, it is natural to suspect that repulsive 
interactions will produce ferromagnetism (FM) due to 
exchange forces, in close analogy to quantum Hall FM 
in flat Landau-levels[19]. Here, we confirm this expecta- 
tion, and show that whereas repulsive interactions disfa- 
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vor d+is pairing, arbitrarily weak repulsion favors FM at 
the edge (even if there is no tendency towards FM in the 
bulk). We identify tunneling-signatures of the edge-FM, 
which can distinguish it from previously proposed d + is 
pairing. 

Flat Edge Bands - To set the stage, we review the 
structure of flat edge bands in the absence of interactions, 
as previously discussed in [6] . We start by considering the 
mean-field Hamiltonian of a superconductor with d x 2_ y 2- 
pairing written as a tight-binding model on the square 
lattice. In the absence of an edge, the Hamiltonian can 
be written in momentum space: 

H d =W*l ^o(k)r 3 + W A (k)n] 



-2t (cos k a + cos kb) — n 
-2A (cos k a — cos kb) 
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Furthermore, we ignore the effect of phase fluctuations, 
which will not be important in what follows, and choose 
Ao to be uniform and real. For convenience, we have 
chosen units of length such that the lattice spacing is 
unity, and shifted the chemical potential such that /i = 
corresponds to half-filling. 

In anticipation of introducing an edge along the (11) di- 
rection, we re- write the Hamiltonian in terms of momenta 
along, k = ka + kb , and perpendicular, k± = ka ~ kb , to the 
edge, where k,k± e [— -^=\. In these coordinates we 
have:"Ho = — 2£fc cos(fcj_) — /x, and Ha = 2A& sin(fcj_) 
where t k — 2t cos k and A^ = 2 A sin k. 

Next, we introduce an edge. Though any direction 
misaligned from (10) will exhibit edge bands, for con- 
creteness we consider an edge along the (11) direction. 
Since k^_ is no longer a good quantum number, one must 
move to a real-space description for the y-direction. For 
each value of k, the Hamiltonian is formally identical to 
an effective Id tight-binding chain with hopping, t k , and 
p-wave pairing Afe. 

It is well known [20] that such Id wires with p-wave 
pairing exhibit zero-energy Majorana end-states so long 
as: \2tk\ < and A(fc) =/= 0. In the present case, there 
will be a Majorana zero-mode for each |fc| < cos -1 
k 0. These edge-bands terminate at bulk nodal points: 
k = 0,±A, where A = cos -1 ^. 

For k between the bulk nodes, the Majorana end-states 
operators with momentum k and z-componcnt of spin a, 
take the form [20, 21]: 

Ika « CfcCT< ^) + i(TS & n ( k ) C -k.-a(y) ( 2 ) 



The states at +k are not independent, but rather are 
related by: 7 fe = -^f_ k u v , where j k = ( 7fe1 - 7*4 ) ■ 



Detailed expressions for the wave-function 4>{y) are 
given in [20] (see also Appendix A) . For our present pur- 
poses, the only important feature of these wave-functions 
is their spatial-extent, and it is sufficient to use the ap- 



proximate form 4>k{y) 



-y/(k 



. The confinement length 



of the edge wave-functions depends on k, diverging as 
6 ~ rei and ~ i±A-fcl near ^ ~ 0, ±A respectively, 

A(l-£ - 2 ). 



^ near k ; 



\k\ sfc ~ |±A-fe| 

and falling to a minimum of 
Here £o = 5^ is the bulk coherence length. 

Instability Towards FM - We incorporate inter- 
actions by an on-site repulsive Hubbard term: Hjj = 
\ J2i n i( n i ~ 1) where is the number of electrons on 
site i. First, we focus only on the sub-space of zero-energy 
edge states to analyze the possible competing instabili- 
ties. From this approach, we identify the dominant ten- 
dency towards FM. Next, we support this picture by nu- 
merically conducting a mean-field analysis that incorpo- 
rates both surface and bulk states. Due to the absence of 
quantum fluctuations about the Ferromagnetic ground- 
state we expect the mean-field description to be sufficient 
at low-temperatures, where thermal fluctuations are not 
important. 

Spontaneous symmetry breaking order can endow the 
flat Majorana edge-bands with dispersion. Focusing only 
on types of order that preserve translation symmetry, the 
generic edge dispersion takes the form: 

H m = Y 7fe [m (k) + m(fc) • cr] 7fe (3) 

k 

where m M € R, and er are spin Pauli matrices. 

The transformation properties of M a b = mo<5 a b+m-er a6 
allow us to identify the physical meaning of the various 
terms in H m . Hermiticity requires M = . In con- 
trast, time reversal acts on M as: T(M) = —M\ indi- 
cating that H m ^ necessarily breaks TRS. Spin rota- 
tions, 7 -> Uj = e- ie - CT 7 with U G SU(2) have the effect 
M — > WMU. Finally, under spatial inversion x — > —x, 
M, transforms acts like: V{M) = -<jyM T a y . These 
considerations show that mo and m break inversion and 
spin-rotation symmetry respectively, allowing us to iden- 
tify mo ^ with the d+is wave pairing, and m/0 with 
edge FM (see Appendix C for an explicit derivation). 
We write the Hubbard U term as: Hjj = 

and focus on 



E r 



(ra tir - n it r) - (n t>r + n^ T ) 
configurations that are uniform along x. Decompos- 
ing the spin-density purely in terms of the edge states: 

Ex ( n T,r - n l,r) « E fe ,J0fc(2/)| 2 7l^7fc + where 
(...) denotes the remaining contributions from bulk 
states gives: 



Hr, 



where V kk > 



J2V kkl {^a^) k {^a^) k , + ... (4) 

k.k' 
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Here we 



see that it is energetically favorable to have ( 7 ^o- z 7 ) 7^ 0, 
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i.e. to magnetically polarize the edge. A similar analysis 
starting by re-writing Hjj in terms of the s-wave pair- 
ing order parameter, H v sa + J2kk> V kk' {lh)k, (7 f 7)fc' 
shows, as expected, that it is energetically costly to form 
s-wave pairing at the edge. 

Further insight into the ferromagnetic edge-instability 
can be gained by developing a Ginzburg-Landau free 
energy for the edge-magnetization. Ignoring, for the 
moment, the bulk states, the effective action for the 
edge-Majorana bands (in imaginary frequency) is: SL = 
\ 12 u ,k Tuj,k (~ iuj ) 7^,fc + Hjj. Introducing the Hubbard- 
Stratonivich field mk ~ J^fc' Vkk< (l^ aZ l) k t to decom- 
pose the quartic term, and integrating out the Majorana 
edge-states, we find (see Appendix D): 



ScS = o E mk ( V 



>kk 



E 



m k \ 



(5) 



kk' 



The unusual, non-analytic \rrik\ term arises from the 
edge-bands' singular density of states in the absence of 
FM order. This term favors m ^ 0, even for arbitrarily 
weak interactions. The saddle point solution 4^ = is: 
m k = J2k' Vkk'- The edge state dispersion is: 



U 



\k\ 



A + \k\ log 



Ifcj 
A + |fc| 



(6) 



The above estimates give: TOfc max ~ Ao. However, this 
is an overestimate for several reasons. First, disorder will 
cxtrinsically broaden the surface-bands, smearing out the 
singular density of states and weakening the instability 
towards FM. Second, edge roughness will further sup- 
press the pairing near the edge, increasing the confine- 
ment lengths £fc, and weakening the residual interactions 
compared to the clean case. Lastly, we have so far ignored 
the mixing between surface and bulk states induced by 
the magnetic order. This is justified only for small split- 
ting (m k < A fc ). 

For larger splitting, the dependence of the confine- 
ment length on the splitting must be self-consistently 



taken into account: 



v±(k) 



limiting the effects 



of interactions for large nik- The bulk states repel the 
surface-bands, keeping them inside the gap, even for ar- 
bitrarily large U. These effects can be accounted for 
by numerically simulating the full bulk plus edge prob- 
lem using Eq. 1 in the presence of an onsite repulsive 



Hu = UJ2i n ii n t -l)~-|Ei miC\a z Ci subject to the 
self-consistency constraint rrii = {c\ a a^ b cn,) (note that we 
allow the mean-field parameters to have arbitrary spa- 
tial variation in the direction perpendicular to the edge) . 
Representative results using reasonable values of t, U, Aq 
are shown in Fig. 1. 

Tunneling Signatures - The edge bands and fer- 
romagnetic splitting can be revealed by tunneling [8- 
10, 13, 14]. At zero temperature, each edge-mode con- 
tributes a Lorentzian peak centered at bias eV = ±mfe, 




FIG. 2. A magnetic field further splits the edge bands, 
and corresponding tunneling peak. Panels a) and b) show 
the edge-spectrum (solid lines) and projected-bulk-spectrum 
(gray shaded regions) for parallel and perpendicular magnetic 
field. (c,d) show the corresponding tunneling conductance 
from the edge states, obtained from a numerical mean-field 
analysis. In black, blue, and red curves correspond to H\\ =0, 
7T, and 14T for (a,c), and to H± =0, lOOmT, and 200mT for 

CM. 



with height and width 7*. s» ^ ps jo\k\- Here, 
70 = 7w(0)|r| 2 is the typical level-broadening due to 
coupling to the metallic lead, where ^(0) is the den- 
sity of states of the metallic lead, and T is the lead- 
superconductor tunneling amplitude. We have calculated 
the detailed tunneling conductance based on a mean-field 
treatment of Hjj and Eq. 1. The results, are in agreement 
with the analytic considerations presented above. 

The tunneling conductance from the edge bands, nor- 
malized to the length of the edge is shown in Fig. lc. For 
weak tunneling, 70 <C 1, the edge-state contribution to 
the tunneling conductance at bias voltage V is given by: 



g(eV) 



2e 2 7o L e k(eV) 
h 2^0 v e (eV) 



(7) 



where L e is the length of the edge, k{e) is the momentum 
of the edge state with splitting e (defined by m&( e ) = s), 
and v e (e) — |^r|m fc =e is the corresponding edge veloc- 
ity. The conductance exhibits a peak near the maximal 
energy splitting, TOfe max , where v e flattens out providing 
large edge-mode density of states that have compara- 
tively large weight near at edge. The width of this peak 
is proportional to 79. 

Parallel Magnetic Field - An applied magnetic field 
H induces further TRS breaking perturbations, and fur- 
ther splits the edge states (see Appendix F for more de- 
tails). We first consider the simpler case of an in-plane 
magnetic field, which affects the edge states through 
the Zeeman coupling: H z s» ^f-H. ■ X)fc7fc <T 7fc- I n the 
FM scenario described above, this Zeeman energy simply 
adds to the spontaneous zero-field splitting, further split- 
ting the edge-bands and tunneling peak (see Fig. 2a, c). 
In fact, such a splitting of the tunneling peak with an 
in-plane field, with slope equal to the Zeeman splitting 
has been observed[12]. 
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Edge Ferromagnetism d+is Edge Pairing 




FIG. 3. Schematic dependence of tunneling peak position, S, 
on magnetic field, for edge-FM (panels a,b) and d + is (panels 
c,d) scenarios. Only the positive bias peak is shown here (neg- 
ative bias is symmetric) . The top (a,c) and bottom (b,d) rows 
show the effects of perpendicular and parallel fields respec- 
tively. The dashed line denotes a first-order phase transition 
where the spontaneous edge-FM is destroyed by the applied 
field. 8z and 8a are respectively the splittings induced by 
the Zeeman energy and orbital screening currents from the 
external field. 

In contrast, no such splitting is expected for the previ- 
ously suggested scenario in which the zero-field splitting 
is due to d+is pairing near the edge[13, 14]. Spontaneous 
d + is edge pairing corresponds to an m type mass (see 
Eq. 3) at zero-field. Applying an in-plane field induces an 
m type term, which further splits the tunneling peak into 
four. However, the further splitting occurs symmetrically 
about the zero-field peak position, and there will be no 
shift in the average peak position. Consequently, the in- 
plane field data of [12] provides evidence for spontaneous 
FM rather than a sub-dominant TR breaking pairing. 

Perpendicular Magnetic Field - In addition to Zee- 
man effects, a magnetic field, H±, perpendicular to the 
ab-plane will produce orbital screening currents. These 
introduce the edge perturbation: Ha = / A • j 
^2 k "io(fc)7fc7fe- Here A = \j,He~ v l Xh x is the vector po- 
tential corresponding to uniform perpendicular field H 
(in the unitary gauge), and m (k) = J dyA(y)\4> k (y)\ 2 w 
jl~A(k). \l is the London penetration depth, H c = 

„£ o \ L is the thermodynamic critical field, and $o = ff 
is the superconducting flux quantum. For YBCO, H c ps 
IT. Note that, in the presence of vortices, A should be 
replaced by A — in the above expression, where 8 

is the superconducting phase. 

For H < H c , the screening currents from H±, induce 
an mo-type term that, in the FM scenario, shifts the 
k > and k < edge-state energies in opposite di- 
rections (see Fig. 2c), leading to a four-fold split tun- 
neling peak (see Fig. 2d). As H is increased beyond 
jrjedge ^ 250mT, the negative energy spin-down states 
are pulled above the chemical potential, and the spon- 
taneous edge FM is killed in a first-order field-induced 
phase transition. Past this point, the tunneling conduc- 
tance shows only a two-fold split tunneling peak, centered 
at eV ~ ±jy-Ao. Simulated tunneling conductance from 



a mean-field treatment of Eq. 1 with Hjj and Ha,z agree 
with the schematic picture presented in Fig. 2. In ad- 
dition, the perpendicular field simulations show a very 
broad peak for H > H° d&e (see Appendix F). 

In contrast, for the previously suggested d + is sce- 
nario, the orbital currents simply add to the zero-field 
splitting. Therefore, in principle, the effects of a per- 
pendicular field on tunneling conductance can provide 
further evidence to distinguish the FM and d + is scenar- 
ios. There is insufficient data from [11, 12] at low fields, 
to discern whether there is first order transition in per- 
pendicular field at H± = H° dsc . It is also possible that 
our estimate of H° dsc is inaccurate due to vortex pin- 
ning effects. There are other complications in the data 
which are not accounted for in our simple model. The 
peak shift saturates at large H±_ , however the previously 
proposed explanation based on non-linear saturation of 
the Meissner currents[14, 22] is also applicable to our 
scenario. There are also complicated hysteresis effects 
in perpendicular field which are not easy to explain [24]. 
For these reasons, we suggest that the parallel field de- 
pendence provide simpler, and more easily interpreted 
evidence. 

Discussion - We have shown that the topological Hat- 
bands provide a useful perspective for discussing TRS 
breaking at the edge of high-Tc cuprate superconduc- 
tors. This viewpoint naturally suggests an instability 
towards ferromagnetic order from repulsive interactions. 
The magnetic-field dependence of tunneling peak split- 
ting provides a simple, though indirect, experimental 
test to distinguish the proposed FM order from previ- 
ously proposed d + is pairing[13, 14]. Alternatively spin- 
polarized planar-tunneling or STM measurements would 
enable one to directly detect of the predicted FM order. 

In closing, we remark that flat Majorana bands are 
also expected to appear in certain classes of nodal spin- 
triplet superconductors [7]. In contrast to the d-wave 
case discussed here, these edge bands would not be spin- 
degenerate, and more complicated types of (incommensu- 
rate density wave) order are required. Moreover, unlike 
the ferromagnetic case discussed here, quantum fluctu- 
ations would generically destroy the density-wave order 
leaving only power-law correlations. 
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Appendix A. Edge State Wave-functions 

As shown in [20], the edge-state zero-mode wave- 
function at fixed k along the edge is: 



A± 



1 



(A+(fc)» + \-(k)v) 



y Nk 

-H± yjy? + 4A(fc) 2 -4t(k) 2 
2(|t(fc)| + |A(fc)|) 



(8) 



where N k — J2 y \^+(^) v + A_(fc) y | 2 ensures proper nor- 
malization. 

A.l - Divergence of Near Bulk Nodes 

Near fc sa 0, A fc -» 0, A± sa V** 2 ^ 2 > • e | A± | = 1. 

Since the wave function decays as e a ln A , this indicates a 
diverging confinement length. Expanding near A m 
gives: 



IA+1 2 



1 - 



2A* 



2(tk + A fc ) ffe 
or ^ = -ln\\ ± \*^*%\k\ i.e.: 



fe^O 



1*1 



(9) 



(10) 



where £o = W^o is the bulk coherence length. 

Expanding near \k\ ~ A, on the other hand, we denote 
tk = n/2 + St k with St > 0, then A + ->• ~^Ak is re S ular 
as <5t — > 0. The diverging binding-length therefore comes 
from |A_| -» 1: 



-A_ 



+ 2VA 2 - fi5t 



/i + 2A + 2Ji 



(11) 



A?( goes to zero rapidly for small y, the wave-function 
can be approximate as |A^| 



y In A_ 



yln(l-5t/A) 



e -y5t/A^ Near fc w ±A; ftfe w 2 isin(fc)c5fc, and A fc w 
2Ao sin(fc). The confinement length in this limit is then: 



£|fc|ssA 



A 



(A-|fc|) 



(12) 



A. 2 - Interpolating Expression 

More generally, 

dE(k,k 



E—Jk) wnere V F ,±(k) = 

k ± =k F ±, and E min (k) = min k± E(k,k±), which 



is the natural coherence length scale for the effective ID 
chain at momentum k. 

Near k w ±A, Vf,± ~ «A = 2Ao, and -E m i n ~ i>f,ii = 
Up = ^/4t 2 — ^ 2 w 2f (for small doping). Conversely, 
near fc « 0, then «f,_l w «f, and Ufji ~ «a- 

Interpolating between these two limits, we expect a 
minimum confinement length of (£fc) m : 



-Jjr near k 



± ^1 — A. Since typically, £o ^ lj (see Fig. 1), there 



are very few states in the region ^1 — < \k\ < 

A. These few states make only a small contribution to 
the energetics of the edge, and for most purposes it is 
sufficient to approximate by the monotonic function: 



Co 



(13) 



Appendix B. Action of Symmetries on 
Edge-Splitting Terms 

Using the convention k > 0, and A real throughout, 
the zero-mode wave-functions can be schematically writ- 
ten as: 



:J 



7fet ~ c kt + ic 



-ki 



7-fet ~ c -fet - ic li 

,t 

fet 



7-fcj. ~ C -H + ic 



(14) 



B.l - Spin Rotation 

Consider a generic (spatially uniform along x) mass term: 



fe>0 



lkaM ab (k)~1- kb 



(15) 



where, a, b e {t, 4-} label spin. Note, we could also extend 
our summation interval to include k < 0, but in this case 
one can easily see (by changing variables k — > — fc and 
then anti-commuting the 7's) that only the components 
satisfying M(—k) = — M T (k) survive the summation; 
consequently the ±fc contributions are not independent. 

Since, "/± ka transform like conventional spinors under 
spin- rot at ion, we see that under 7 — »■ e l8 ' rT ^f = Ue-f, we 
have: M -> Uj MU e = &»\J- x a y MU. This transforma- 
tion suggests that we parameterize A'/ = a v M 1 such that 
M transforms as a (1/2, 1/2) tensor under spin-rotations. 

B.2 - Hermitian Conjugation 

Under Hermitian conjugation the edge states transform 
as: 



7fe = °lb1-kh 



(16) 



Correspondingly, the edge-splitting terms transform like: 

(lkaMab(kh- kb y = K b (kh-kb^ka 

= K b (k) (-^7 fe6 0K y a<7- fea <) 



Ika 



ab 



7-feb 



(17) 



This constrains M — M\ which can be accomplished by 
parameterizing M = tuq + m • er. 
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B.3 - Time-Reversal 

Time-reversal (TR) acts on the electron operators like 
Cfea — > i&yKckb where K represents complex conjugation. 
In terms of the edge modes, one then has ^ka —> i^ab^-kb, 
and: 

(lkaM ab (k)-f-kb) -> i(T v aa n-ka'M* ab i(T v bb rtkb' 

= lkb>{o v )l b Miy aall _ ka , (18) 

or equivalently: 

M -M f (19) 

We see that the constraint of Hermiticity necessarily re- 
quires that any mass term break TR symmetry. 

B.4 - Spatial Inversion 

Inverting the x-direction acts projectively in the presence 
of d-wavc pairing. In order to compensate for the sign 
change of the order parameter under x — > —x, we must 
accompany this inversion by a global gauge transforma- 
tion c — > ±ic. Since the particular sign is irrelevant, let 
us choose c — > ic for definitencss, in which case we have 
7 ±fecr — > ij T kcr- Note that, in contrast, d + is pairing 
transforms to d — is pairing under this transformation. 
It is not possible to relate d ± is pairing by a global 
gauge transformation as the relative sign of the d- and s- 
wave components is physically observable, corresponding 
to the direction of edge currents from the TRS breaking 
pairing. 

The mass term then transforms as: 

M -> -a v M T a y (20) 

Using the parameterization M = m + m-er, m — > — ttiqI 
and m > m. 



Appendix C. Explicit Representation of Electron 
Bilinears by Edge States 

Decomposing the electron plane-wave operators Ck a 
into the single particle eigenstates we have: 

c k\(y) = -iMvh-u + ■■■ 

cfet(y) = <Afe(y)7fet + ••■ 

c li(y) = i<t>k(vh-kt + ■ ■ ■ 

cki(y) = 4>k{y)iki + • • • (21) 

where we have explicitly written only the contribution 
from the zero-mode edge states, and (...) indicate con- 
tributions from other extended states, which we ignore 
for this section. Similarly we may decompose the spin- 
density and s-wave pairing in terms of single-particle 



eigenstates. Retaining only the contributions from the 
zero-energy edge states we have: 

S z (q = 0) = E [4f c fcT - C /4 C *4 

k 

~ E l^' 2 H7-/47fcT _ (h-kflki)] + ■■■ 

k 

= E l^ fc | 2 i^ktl-ki + «7fe4.7-fef] + • • • 

k 

= Ei^N ct ^ + --- ( 22 ) 

k 

Similar, the s-wave pairing density is: 

T ' = E ( C lt C -fe4- + C -kl c kt) 

k 

~ E ^k [(-*7-fc;) (*7fct) + ft-ki) (7fcf)] + ■ ■ ■ 

k 

= -E^Hl7fe + --- (23) 
k 

Note that, for repulsive on-site interactions Urn{rii — 
1) can be decomposed as —US\S l z or U\T S \ 2 , indicating 
that the m terms is energetically favored by repulsive 
interactions (U > 0), whereas the mo mass is not. 

Appendix D. Derivation of Effective 
Landau-Ginzburg Action 

Consider just the edge modes, then the effective action 

is: 

S = Y 7^7 - \ E Vkk' ( 7 V 7 ) fc (M), (24) 

Applying Hubbard-Stratonovich to decouple U term: 

5 = 2 S m ^fcfe' m 'fe + E T f + mfcfT ") 7 (25) 

Integrating out 7's gives a term: 

5 7 = trln [iuj — mu<j z ] = — ^ In [cj 2 + m 2 ,] + const 

uj,k 

(26) 

The w integration can be performed by subtracting an 
m independent constant to regulate the uj — >• oo behav- 
ior: J2 U ln (^ 2 + m fe) XL ln ^ Here > arbitrary 
(for concreteness, we take S < m), and we are only in- 
terested in the m dependence of the final answer (since 
m-independent terms will not affect the saddle-point so- 
lution). The term ln (w 2 + m 2 ) has a branch cut, which 
can be taken between w = ix, x G [— |m|, |m|]. Simi- 
larly, the ln (w 2 + S 2 ) branch cut can be taken between 
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[—iS, iS\. The branch cut of the two logs cancel between 

2. 2 

co E [—i5,iS] when taken together as In ^2 + ™2° , leaving 
the integrand well defined on the real axis. Moreover, 
since as \oj\ — > oo, the integrand falls off as ~ we 
may extend the integration range to a contour encir- 
cling the full upper half-plane. Since, in the upper half- 
plane (UHP), the integrand is analytic everywhere out- 
side bj G [iS, im] on the imaginary axis, the contour can 
be shrunk to one encircling (but infinitesimally outside) 
the branch-cut in the UHP. On each side of the branch 
cut, the integrand differs by — 2iri (for u — xe l ^^ 2±e \ 
with x G [S, to], then as e — >• In (w 2 + to 2 ) is in- 
dependent of e but — In (cj 2 + S 2 ) differs by — 2ni be- 
tween +e and — e). Consequently the integral reduces to 
— 2tti f-g™ 1 ^ dui = \m\ — S Consequently we find: 

5 7 = - \m k \ + const (27) 

k 

A less direct, but much simpler way to evaluate the 
sum is to note that that J^- = J2u> J+m 2 = s g n ( m *0' 
which indicates that S 7 = — |mfe| + const. 

For repulsive interactions (i.e. Vkk 1 > 0), we can choose 
to > 0, and drop the absolute value signs. Therefore, we 
see that in the free energy for to, there is a linear term 
due to the singular density of states of the flat bands in 
the to — y limit. The resulting saddle point solution and 
energy are simply: 

m k |s P = ^2 Vkk' S\ sp = - ^ 51 Vkk ' ( 28 ) 

k' kk' 

This consideration is useful, because, generically the 
bulk will add a term that suppresses mk, which we 
can write schematically as J2 r \ rn k\ 2 - However, due to 
the linear term from the edge states, the saddle point 
will always have non-zero to, and energy that tends 

as ~ kk r2 kk ' even when magnetization is arbitrarily 

strongly disfavored in the bulk (r ^ oo). 

D.l - Mean-Field Energy and Dispersion 

The end-state wave-functions may be reasonably approx- 
imated by: <f>k(y) ~ -y=e~ v ^ k . This ignores the oscilla- 

V sfc 

tory structure, which is reasonable when calculating ma- 
trix elements like / dy\<ftk{y)\ 2 \4'k'(y)\ 2 that will appear 
below, since, for k very different from k' , the oscillations 
in (f>k and <pk' happen with very different periods, in which 
case \(j)k\ 2 sees roughly the average of |0fe'| 2 . 

At zero-temperature, approximating w £o/|&| ; the 
The dispersion m k is approximately given by: 

= ^ k ( A+klo ^khd) (29) 



This is strictly only valid near k w 0, and that as k 
becomes very close to ±A, the spectrum should dip back 
to zero. To capture this behavior, we should replace: 
k — > ^ 2 (A — |fc|) in the above expression. However, for 
energetic purposes, most of the states have dispersion 
following that described above. 

Similarly, the mean-field energy gain is: 
£mf = E »W = ^ jf dkk (a + * log (^) ) 



Appendix E. Details of Mean-Field Simulations 

Mean-field simulations of the tight-binding model Eq. 1 
of the main text, with an on-site Hubbard U term were 
conducted numerically. We decomposed the Hubbard U 
interaction in terms of spin, and considered magnetiza- 
tion profiles that uniform along the edge (x-direction) , 
but spatially inhomogeneous along y. In this case, trans- 
lational symmetry was preserved along x, and the states 
at different k were decoupled, enabling the problem to be 
reduced to a collection of separate ID models with the 
same magnetization profile. The simulations displayed in 
the main text discretized k into 51 values (the results did 
not change markedly for larger numbers of k), and the 
system length in the y-direction was L y — 900 sites (fur- 
ther increasing the length led only to minimal changes). 
Self-consistency was achieved by starting with an initial 
ansatz for the magnetization, and iteratively computing 
the induced mean-field magnetization until the fractional 
(root-mean-squared) change in magnetization between it- 
erations was < 1%. 

Tunneling conductance was computed as ^-tr|r e /,,| 2 , 
where r e h is the reflection matrix for processes in which 
an incoming electron from the lead is Andreev reflected 
as a hole. We assumed a uniform coupling to between the 
lead and the wire, that is independent of the momentum 
along the edge. This is appropriate for a thin, high tun- 
neling barrier. Since the effective barrier thickness de- 
pends on the incident angle of the incoming electron, for 
a thicker barrier, the lead-superconductor coupling will 
depend on the incident angle of the incoming electron. 
This will obscure the underlying behavior of the tunnel- 
ing peak by providing a strongly momentum dependent 
visibility to the edge-states. 

In the next appendix, we provide further details about 
the impact of perpendicular magnetic fields, and incor- 
porate orbital magnetic field effects into the numerical 
mean-field simulations. 
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Appendix F. Orbital Magnetic Field Effects 

This appendix further details on the orbital effects of 
a perpendicular field. For H < H c \ = w 30mT, out 
of plane fields are screened by the superconductor as: 



H(y)z = Hze- y ' XL 



(31) 



where Xl is the London penetration depth (Xl 3> £ 
for cuprates). For H > H c \, vortices enter the system 
and the current profile at the edge depends in detail on 
the precise distribution of vortices, which will depend on 
pinning effects from the sample edge and from impurities. 
A realistic description of perpendicular field effects in 
actual cuprate edge-tunneling experiments would need 
to confront this complication. However, it is useful to 
gain some intuition from the simpler case of uniform edge 
currents, and we will imagine throughout that the field 
is below H c i (even though we consider fields that would 
exceed the H c \ for realistic cuprate materials). 

Working in the unitary gauge (where the supercon- 
ducting phase is fixed uniformly to zero), H can be pa- 
rameterized by the corresponding vector potential: 



A = HX L e- y/XL x 



(32) 



where Xl is the bulk penetration depth. This induces 
the perturbation (that we re-write in terms of the edge 
states): 

J A 'J = e X! J d V A x{y)v x {k)c\ v c ky 



eBX L v F ^irlhk 



itBXlZq 



B 



J2 A ( fe hhfc 



(33) 



where in the first line, we have used £ -C Xl, which 
means that for most k, <C (except those very close 
to the bulk nodes, which contribute only weakly to tun- 
neling conductance). In the last line, <3>o = ft = f * s the 
superconducting flux quantum, and we have introduced 



the critical field H c = 



$0 



(w IT for YBCO). 



F.l - Orbital Currents in the Vortex State 

In this section, we consider the effect of a finite density 
of vortices on the edge-current profile. Generically, we 
find that the typical current density near the edge is un- 
affected by the presence of vortices. 

In the extreme type-II limit, valid for the quasi-2D 
high-Tc materials, the vortex cores do not strongly effect 
the current distribution at the edge. Neglecting the con- 
tribution from vortex cores, the Landau-Ginzburg free 



energy density of the superconductor in a perpendicular 
external field, H is: 



Am 



V6> 



2e 



(V x A - Hz) 2 (34) 



where p s is the superfluid density, 9 is the phase of the su- 
perconducting order parameter, and the Hz term comes 
from external current sources that attempt to impose a 
constant magnetic field: B = H. 

It is convenient to absorb the phase-gradient, V6> into 
the vector potential: A — > A = A + j-V6. It is then A 
which acts as the orbital perturbation for the edge states 
rather than A. So long as we impose A • h = 0, where 
h is the normal vector for the sample boundary, then we 
can interpret A as being proportional to the supercurrent 
j s = — |^A. It is crucial to recognize that it is A and not 
A which couples to the edge states. We shall see that, 
while A is strongly dependent on the vortex distriubtion, 
near the sample edge, A is not. 

With this transformation the free-energy density reads: 



4tt^= J_|A|: 



+ - V x A + $oPtj(r) - H 



(35) 



he 



IS 



z-vx(ve) 



where Xl is the London penetration depth, $ 
the superconducting flux quantum and p v (r) — 
is the density of vortices. Choosing A = A x (y)x, and 
minimizing with respect to A, gives the equation of mo- 
tion: 



-d, 



1 



A x {y) = -§ d y p v {y) 



(36) 



We approximate the vortex distribution by a uniform 
continuous density (valid in the limit where the inter- 
vortex spacing is much shorter than the Xl, but much 
larger than the core size £o)- Specifically, we consider 
p v (r) = -^9{y—L) where 9 is the step-function, and L is 
some unknown distance, which can, for example, account 
for barrier effects repelling the vortices from the sample 
edge. 

Solving the above equations, subject to the boundary 
condition that —d y A x (y = 0) = H gives: 



A = HX L x [ e -lf- L l/^ + (1 - e - L ' XL )e 



(37) 



Since the edge states reside predominately within dis- 
tance < ^ of the edge, they are effected by A x (y — > 
0) = HXl- This result is independent of L, indicating 
that the typical size of the orbital perturbation is insen- 
sitive to the presence and distribution of vortices. 

Note that, the current profile will be more complicated 
for the more realistic case of a discrete vortex lattice, 
pinned locally disorder. However, by the above consid- 
erations, we expect the overall magnitude of the super- 
currents at the edge to be roughly independent of these 
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FIG. 4. Tunneling peak splitting in a perpendicular magnetic 
field for L y = 900, t = 240meV, A = lOmeV , U = 2eV, and 
7o = 0.4Ao. Zeeman splitting is neglected and only positive 
bias is shown for simplicity. At zero field, the peak is split 
by the spontaneous edge FM. For H < Hf sc « 0.26T (top 
panel), the zero field peak splits further due to the screening 
currents induced by the orbital field. The spontaneous FM is 
killed in a field induced first order transition occurs at H = 
Hf £C . For H > Hf sc there is only a single split peak. 




0.2 f 0.4 0.6 0.8 
H ed g e H ^ (X) 



FIG. 5. Color plot of the field dependence of tunneling con- 
ductance, corresponding to Fig. 4. 

details. In this light, it is difficult to understand how 
the presence or absence of vortices can account for the 
substantial saturation and hysteresis effects observed in 
[12]. 

F.2 - Tunneling Conductance in Perpendicular Field 

Simulations of the tunneling conductance, due to An- 
dreev reflection from the edge-states, in the presence 
of a perpendicular field are shown in Figures 4 and 5, 
for the case where there is spontaneous edge FM at 
zero-field. The top panel of Fig. 4 shows results for 
H < Hfs c w 260mT, where the initially FM-split peak 
is further split into two by the supercurrents induced by 
the external field. For larger fields, H > iJ° dgo , the zero- 
field FM splitting is discontinuously destroyed by the ap- 
plied field, resulting in a single split peak (see bottom 
panel of Fig. 4). At high-fields, the single peak continues 
to move out with increasing H±, and broadens substan- 
tially. This is in reasonable agreement with experiments, 
and the substantial broadening may obscure the visibility 
of the predicted low-field first order transition. 

The corresponding tunneling data is shown in a false 
color plot in Fig. 5. 



